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General Abstract 

This documentation presents an input-output model which has been modi- 
fied to include the environmental impact of economic operation. In lieu of 
market prices for the environmental factors, trade-offs with regional income 
and employment are estimated for use in regional planning. The program is 
written in FORTRAN IV with single precision for the IBM 360/65 system. The 
example data has been set-up to contain five endogenous sectors and three 
exogenous sectors (including households). The data is then compressed to 
three endogenous sectors. 
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PURPOSE AND GENERAL CAPABILITIES OF THE PROGRAM 



Input-Output is a well known analytical tool which is particularly well 
suited for gaining information on the economic aspects of environmental 
quality and resource use. These economic aspects concern trade-offs between 
levels of pollution or resource use and levels of output, income, and employ- 
ment in the economy • When an environmental sector is added to the I-O matrix, 
these trade-offs may be estimated by using interdependence coefficients and 
the multiplier concept normally incorporated in input-output models. The 
purpose of this documentation is to provide information about a computer pro- 
gram (ECON-ECOL) that allows an analysis of economic and environmental data 
via tine input-output approach by estimating interrelationships and multipliers 
which link the economy with its demand for resources and its supply of pollu- 
tants. , 

The details of the program's operation are described in later sections, 
but a statement of the general aspects of the program is useful as an intro- 
duction. The program links the environmental factors with the economic sec- 
tors by expressing resources as inputs per dollar of- output and pollutants 
as production by-products per dollar of output. Environmental factors are 
introduced into the input-output model as L sectors of positive environmental 
imports, (resource inputs) and K sectors of negative environmental imports 
(exported pollutants) . The result is a MXN (where M = K + L) environmental 
matrix. The data in this matrix are converted to quantities per unit of 
output (production coefficients) for each economic sector. This yields a 
more complete production funct-^'^n for each economic sector of the model. 

Once formed, the env^ coefficient matrix is post multiplied 



by the inverse of the economic processing matrix. The resulting matrix con- 
tains environmental-output (E-0) multipliers that indicate the direct and 
indirect effects of economic production on the environment. By using the 
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inverse of a processing sector containing households, the induced environ- 
mental effect is incorporated in the E-0 multipliers. Finally, by combining 
the elements of the E-0 matrix with appropriate income and employment data, 
environmental trade-offs in terms of dollars of income and/or levels of em- 
ployment are generated. 

The program is written to provide maximum information and flexibility 
and to be relatively inexpensive with respect to computer operation. The 
information provided as output includes multipliers calculated with house- 
holds both exogenous (type I) and endogenous (type II). The flexibility 
provides for any size matrix up to 95 x 95 for the complete I-G model, up 
to 90 x 90 for the processing sector, and up to 25 environmental factors. 
The environmental factors can be either pollutants or resources or both. A 
compression subroutine allows a large model to be aggregated to as small a 
matrix as is desired. Other options allow for deleting the emplojmient vec- 
tors or for operating simply as an I-O model excluding the environmental analys 



USER DOCUMENTATION 

Th^t use of the various options of the ECON-ECQL program results in 
differing data requirements, arrangements, and coding. The use of the option 
COMPRESS calls for a iilight alteration in the way in which information is 
presented to the computer. Information on data requirements and the arrange- 
ment of data relative to the main program is presented here as if the pro- 
gram '^s full capabilities are to be used. 

Data Requirements 

Assuming that an analysis including employinent-environmental relation- 
ships is desired, data is needed for two matrices and two vectors* The 
matrices are for economic data on each industry's sales and purchases and 
environmental data on resources* required and pollutants produced by each 
industry listed in the processing sector of the I-O model. The two vectors 
consist of employment totals for each industry and a list of names for the 
industries and environmental factors that constitute the production function - 

Economic Data 

It will be assumed that the user is familiar with the data requirements 
of the economic sector of this input-output model. Therefore, the only dis- 
cussion of this topic center on level of aggregation for each sector. 
This topic will be taken up later. 

Environmental Data 

The environmental data must contain environmental inputs for each in- 
dustry of the processing sector expressed as positive quantities, and 



pollutants expressed as negative quantities. The number and type of these 
positive and negative environmental factors depends on existing clef initions ^ 
of pollutants, data availability, and the user^s analytical interests. Tlie . 
environmental data should conform to the economic data in terms of indus- 
trial classification, production period, etc. However, unlike the economic 
data that is expressed in dollars, resources and/or pollutants may be express- 
ed in any appropriate units (e.g. tons, gallons, acre feet). Differing units 
may be used within the same environmental matrix. 

Employment Data 

Employment totals for each of the industries within the processing sec- 
tor are used to calculate environmental-employment trade-offs as well as the 
standard I-O employment multipliers. 

Row Names . 

For print-oUt purposes row names for the economic and environmental 
matrices are needed. The economic row names include processing industries, 
households, and final purchases. Names for final demand sectors are not 
needed as they are listed as columns which are numbered rather than named. 
Environmental row names contain designations for each resource and each 
pollutant for which data has been collected. 

Each of the above data items may be thought of as separate data decks 
in the program. Hereafter, they are referred to as ECQN for the economic 
data, RES (resources) for the positive and negative environmental data, 
EMPLOY for employment, and NMi E for row name vector, 

(Control Ca rds 

'L'lic final data category contains control and option cards. The options 
desired arc listed on one computer card, called the option card (designated 
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OPT). The options available are for environmental analysis, employment- 
.environmental analysis, and compression of original data. If the environ- 
mental analysis is excluded, resource and pollutant data is not needed, and 
the prograiTi will perform a standard Input-Output analysis. This analysis 
results in output multipliers, type' I & II jLncome multipliers, and if de- 
sired, employment multipliers. Employment multipliers for the economic 
sectors and the environmental sectors can be excluded via the second option. 
Through the compression option it is possible to reduce the matrix con- 
structed from the original data collection to as small a size as is desir- 
able. The result will be a reduced economic and resource matrix with a 
correspondingly shortened employment vector. 

The control card SIZE specifies the size and degree of partitioning of 
the matrix that is to be used in the analysis for that particular computer 
run. The size of the matrix refers to the entire transaction table. The 
partitioning consists of separating and finding subtotals for the process- 
ing matrix, the value added vectors, and the import vectors. 

The final set of control cards are used only if subroutine compress is 
called. These cards describe which rows and columns are to be aggregated 
and the position of the resulting rows and columns in the reduced matrix. 
These control cards, then, constitute another data 'Meek*' and will be desig- 
nated MOVE . 

Generalized Job Stream 

The control card, option card, and data decks are combined with the 
main program as illustrated in Figure 2. 
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Data Coding 



Industrial Classification 

The level of aggregation for ':hich the data is collected depends on the 
industry classification, Thiii classification also determines the size of the 
matrices to be used, ihis program is written in general terms so that any 
size matrices can be used that are less than or equal to the maximum speci- 
fied in Ltie dimension statements. These limits consist of 95 x 9.^ for the 
JO transaction table, 89 x 89 for the processing sector without households, 
and 25 x 100 for the resource matrix. Larger or smaller limits can be used 
but it will be necessary to change several cards in the program. This infor- 



mation can be found in the section describing the technical details of the 
program itself. 

Economic and Resource Format 

The collected data must be coded so that they can be read by the compu- 
ter. The economic and environmental dar.^ are presented in the same form. 
Each quantity to be placed in a matrix is accompanied by the row and column 
number which designates its position in the matrix. Each data card will 
have four data points. Each data point will consist of three numbers. These 
numbers are row, column, and data value. The exact format is specif ed as 



4(213, F12.3), 



Card Item 



C^rd Column 



Format 



M 1) Row number 
2) Column number 



4-6 



1-3 



13 



13 



3) Datum value 

4) Row number 



19-21 



7-18 



F12,3 



13 



5) Column number 



22-24 



13 



V 



9 



Card Item 



Card Column 



Format 



6) Datum value 



25-36 



F12.3 



10) Row number 

11) Column number 

12) Datum value 



55-57 



61-72 



58-60 



F12,3 



13 



13 



This arrangement of data is repeated on as many additional cards as needed. 
Data Arrangement 

Zeroes need not be coded since all matrices are filled with zeroes ini- 
tially and whatever cells are not filled with data will automatically re- 



the data within the decks ECON & RES is unimportant. Each deck could be 
shuffled and the row and column designation will insure proper placement of 
the quantity. 

There are only two restrictions on the arrangement of the data. First, 
the households sector must be the vector Immediately following the process- 
ing matrix* This facilitates bringing households inside the processing sec- 
tor for estimating type II multipliers. It is, therefore, necessary for the 
household row and column total to be balanced. Second, all value added sectors 
must be grouped together and preceed a similar grouping of the import sectors. 
These groupings are necessary because of the feature that partitions the 
transaction table and presents row sub-totals for each of the groups mentioned. 

Pollutants are treated as negative resources. Therefore, each pollu- 
tant quantity should be punched as a negative number. The last card in deck 
ECON and deck RES must be -99 starting in column one. This tells the compu- 
ter that the last card in that deck has been read. 

Employment Format 

Employment totals are placed six to a card. Each figure is right 



gister as zero. An additional aspect of this feature is that the order of 
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justified in a ten space field beginning in column one* As many cards as 
are necessary may be used* 



Item 

Card 1 

1. Employment total sector 1 

2. Employment total sector 2 

6. Employment total sector 6 
Card 2 

1. Employment total sector 7 

2, Employment total sector 8 

6. Employment total sector 12 



Card Column 

1-10 
11-20 

51-60 



Format 
110 

no 

110 



1-10 

11-20 
» 

51-60 



110 
110 

no 



Row Name Format 

Row names are punched one to a card. The first three spaces contain the 
row number. The next three spaces are left blank, and the names are then 
placed anywhere within a twelve character field- This is done for both eco- 
nomic sectors and environmental categories. 



Item 

Card I 

1. Row number 

2. Row name 

Card 2 

1- Row number 
2. Row name 



Card Column 

1-3 
7-18 

1-3 
7-18 



Format 

13 

3A4 

13 

3A4 
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Control Card Format, SIZE 

Control uards SIZE and OPT have the appropriate numbers right justified 
in fields of five spaces each. The first card, SIZE , contains four numbers 
relating to the economic matrix being used for analysis. These numbers are 
the size of the transaction table (NN) , the number of the last row of the 
processing sector (NP) , the number of the last row of the value added sector 
(NV) 5 and the number of the last row of the import sector (NI) • 



Item Card Column Format 

1, NN 1-5 15 

1. NP 6-10 15 , 

3, NV 11-15 15 

4, NI 16-20 15 



Control Card Format, OPT 

The option card contains three numbers. The first (MN) is for the 
environmental analysis. If environmental multipliers are being calculated 
MN =» the number of environmental categories (positive & negative) that are 
being used. If no environmental analysis is desired MN = 0. The second 
number (LAB) indicates the length of the original (i.e., noncompressed) 
employment vector. If no employment data is present LAB = 0. The last 
number (NCOMP) indicates whether the matrix to be analyzed is a compressed 
version of the data being read in. If compression is desired, NCOMP - the 
number of times that several sectors are being compressed into one. If sub- 
routine compress is not required, then NCOMP = 0. 



Option Card Forma t 



Item 



Card Column 



Format 



1. 



Environmental Option (MN) 

MN = 0 deletes option 

M]^ = number of environmental 



1-5 



15 



sectors calls option 



2. 



Employment option (LAB) 
LAB 0 deletes option 
LAB = number of processing 



6-10 



15 



sectors calls option 



3. 



Compression Option (NCOMP) 
NCOM? = 0 deletes option 
NCOMP = number of compressions 



11-15 



15 



within matrix calls 
option 



If NCOMP ^ 0, the size of the matrix being used for analysis is smaller than 
the matrix being read in. Therefore, the data on card SIZE (NN, NP, NV, NI) 
refer to the compressed matrix and are smaller than they would be if compres- 
sion were not desired. 

Control Card Format, MOVE 

The last set of data cards (deck MOVE) provide the information needed 
by the compression subroutine. The information consists of the sizri of the 
transaction matrix being read (since compress is being called, this number is 
larger than its counterpart on control card SIZE ) and a vector of numbers in 
groups of three which states the destination, origin, and range of each com- 
pression. 

Only adjacent sectors can be aggregated, and while some column and row num- 
bers will change, the sequence must be maintained. The first group of sectors 
which are aggregated are placed in the row and column which the first sector 
of that group occupied. The unaggregated sectors between the first and second 
compression will automatically be moved over, and the sum of the second group 
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of aggregated sectors will be placed in the following row and column. The 
information presented in the compression vector, then, is (1) the number of 
the new column (row) which contains the aggregated sectors of the initial 
compression, (2) the number of the first sector in the group to be added 
(i.e. , the starting position of the compression) and (3) the number of the 
last sector of the group (i.e., the end of that particular compression). 
Similar information is given for each compression. One compression in-- 
eludes summation of both columns and rows. 

The data for subroutine COMPRESS are right justified in adjacent fields, 
containing five columns each. The first field begins in column one. The 
first card contains the size of the transaction matrix being read in and as 
many as four groups of compression data. If more than four compressions are 
desired the information is coded cn subsequent cards as needed. All addi- 
tional cards contain only the four groups of three numbers. 

Move Format . 

C ard Item Card Column Format 

Card 1 

1. Dimension of matrix read in 1-5 15 
Dimension = N for an NxN 

matrix 

2. Column (Row) number of 6-10 15 
position of aggregate vector 

3. First column (Row) of group to 11-15 15 
be aggregated 

4. Last column (Row) of group to 16-20 15 
be aggregated 

5. Column (Row) number of 21-25 '15 
position of aggregate vector 



14 



Move Format (continued) 

Card Item Card Column Format 

6. First column (Row of group to 26-30 1-5 
be aggregated 

• • • . . 

• • • 

13. Last column (Row) of group to 61-65 1-5 
be aggregated 

Card 2 

1. Column (Row) number of 1-5 1-5 
position of aggregated vector 

2. First column (Row) of group to 6-10 1-5 
be aggregated 

When COMPRESS is used the row names in NAME correspond to the aggregated 
sectors rather than the original sectors. 

Example Code Sheets 
Examples of the coding, procedures are presented in the following code 
sheets. The data used are from a problem run as an illustration of the 
program's use and output • 
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Complete Job Stream 

Control cards for the computer system consist of a job card, a class 
card, and supporting JCL cards. At maximum a computer run with this program 
should not exceed 3 minutes, 5000 lines, or require more thc'n 320K for storage 
space. It is possible that for smaller data sets smaller parameters may be 
desired in order to achieve faster turn around. 

The program is written for Fortran G. Due to core requirements it 
cannot be run on the special superfast compilers such as WATFIV. JCL cards 
used with the program are indicated in the complete job stream (figure 3). 

Output 

The results of this program contain matrices presenting standard results 
of input-output analysis (i.e., transaction matrix, direct coefficients, and 
inverses with and without households), vectors of income and employment mul- 
tipliers, and vectors of environmental multipliers. 

The environmental multipliers present the trade-offs between the en- 
vironmental . factors and output, income, and employment. In addition, resource- 
resource multipliers are presented which show the total amount of resource 
use in the region for a unit increase in its usage by each industry* 

With the exception of transaction table all output is self-explanatory- 
A possible source of confusion is the numbered columns of this matrix. Within 
the processing section the columns correspond to the named rows. However, 
since final demand categories can be different from final payments, the column 
designations must be kept track of separately by the analyst. Finally the 
last two columns contain respectively sub-totals of the processing sectors and 
the total of the entire row. 
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Example Problem 

The following is output from a sample set of data. The data is from an 
8x8 transactio n matrix whichr jias been reduced ^ via COMPRESS ^ to 6 x 6 . 
It should be noted that it is impossible to get output of data before and 
after compression on the same computer run (i.e., NCOMP cannot = 0 and 0 
at the same time) . What follows is the combination of two runs for example 
purposes. 

The environmental factors consists of two pollutants and one resource. 
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ORIGINAL ECON DATA 



.1 . AG'^ITMLTItpf 

nn(jc,rMni Q<; 
T Monp r s 



I )0.0 00 

0.0 
> "10.000 
"■■OO.OOO 
■■'00.000 
^ )0, 000 
7OC.00O 

?oo,'ooo 



5 0. 000 
7S.00O 

■'OOo 000 

0.0 
^00, 000 
■=^00.. 000 
^ 7'^.»000 
^00,000 



0.0 
^P'l. 000. 
100, 000 
300.000 
^^00.000 
? 000 
lOOo 00 0 

0.0 



•^00.000 
^00.000 

0.0 
.''SO. ^00 
=100.000 
^iSO.OOO 

ion, ono 

TOO. 000 



1, AGRICUL'^UHF 
?. MTNPS 
•^.rONST-MANilF 
^, UTILITT^*: 

cnMMcpfr 

6. MOUS*^HnLD<: 

7, r,OVFOMMFNT 



?. 0.0, 000 
3 ,)0,0 00 

«oo.ooo 

^lOO.OOO 
'\ )0. 000 
^ )0.0 00 
0,0 ' 

0.0 . 



"^00, 000 
750, 000 
100.000 
A 00. 000 
-i 00. 000 

0.0 

0.0 
■^00. 000 



500,000 
0.0 

0.0 I 

0.0 ' 
^00,000 
^P^i.OOO 

0.0 

0* 1 



SSO.OOO. 
750. 000 
0.0 

"«'=;-o.ooo 

0,'> . 
0, 0 
0.0 
0. n 



COMPRESSED OUTPUT 



21 



1, AG» !-MIMr s 
3,irTL-Cnvw 



"5, 00 

^'OO.- 00 
^^0O. 00 



?• 

?.?H,00 
1 00, 00 
7 00. 0{3 



1?0 0.00, 

flOOoon 

l?50.00 



• 4 
1. ??0.00 
00« 00 
ROO. 00 



142^.00 10.?5.00 ?'^50.00 .?!T'^0,00 



1 300o 00 

^ :)75.oo 



PTS.OO 
! 00, 00 



650, 00 
100. 00 



0.0 

- n.n 



■' 375, 00 



375<.00 



750.00 



0. 0 



A, TMnoC;T.< 



. TOTAL TMoocTc 



5 00, 00 



5 00, 00 



0. 0 



0, 0 



•^00,00 



■^00.00 



■»00. 00 

^oo;oo 



' CHLUMN TOTAL <: 



► 100. Oa W'OO.OO 4400.00 . ?450.00 



5 . TRANSACT TPN T^PLT (continued) 



?o ^•"NSt-M^^;|)F 



5 

500. 00 
0. 0 

4 00. 00 



6 

900. 00 
0. 0 

^5 0. 00 



(Sub- total) (Total) 



1 h 50^00 
13 00.00 

.?fl 50. 00 



43OO0 00 
1400.00 
4400.00 



INT. pi)PCHASP5 



500.00 1750.00 5P00,no 10500.00 



4o HOUSPH^LOS 



^:5, 00 
0. 0 



0.0 

0. 0 



.2225,00 
1275, 00.. 



2450.00 
1 275,00 



? ?5, 00 



0. 0 



3500. 00 



BT?*?, 00 



6. IMPORT S 



THTAL IMopfTc 



0. 0 



0. 0 



0, 0 



0. 0 



800.00 



POO. 00 



1 1 00, 00 
1 1 00, 00 



COLUMN TOTM.c 



1125.00 1250. 00 



[ERJCbqw su^ - 140:"^.-) 000 



0,0 



COLUMN ';UM = 



o„o 



14O25.0000 



ENVIRONMEK'TAL FACTOR'^ 



1» SULFIDES 
2. QOD 

i. WATER RFC 



1 

•25C0 00 
-3C0* 00 
3300.00 



2 

0 0 
- 100. 00 

1 100. on 



-flOO.OO 
-^25, 00 
5 00 . 00 



4 

-200.00 
0,0 

200-./ 00 



EMPLOYMtNT VFCTOk 



1. 
2. 

3* 



6 75 , 00 
400.00 
325.00 



DIRECT REeUIPEMENTS (Econ technical coefficients) 



1. AGRI-MIMES 

2. C0NST-MANUF 

3. UT IL-COMM 

4. HOUSEHOLDS 

5. GOVERNMENT 
6. IMPORTS 



1 

0.05233 
0. 09302 
0. 18605 
0, 302 33 
0. 25000 
0. 11628 



0. 16071 
0. 07143 
0. 50000 
0. 196 43 
0, 071 43 
0.0 



0. 27273 
0-, 18182 
0 30682 
0. 14773 
0. 02273 
0. 0681 8 



COLUMN TOTALS 



l-OOOOO 



1 00000 



1 .00000 



INTERDEPENDENCE TABLE 

I 2 3 

1. AGRI-MINFS 1. 23319 0. 55276 0.63018 

2. C0NST-MANUF 0. 21932 1. 35235 0,44100 

3. UTIL-C0MM 0. 48918 1. 12382 1. 929tt6 



COLUMN TOTALS 



1. 94169 



3- 02 893 



3. 00104 
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DIRECT R60UIREMENTS (Econ technical coefficients) 



ImAGRI-MINES 
2«C0NST-MANUF 

3. UTlL-CaMM 

4. HOUSEHOLDS 

5. GOVERNMENT 
6". IMPORTS 



0.05233 
0. 09302 
0. 18605 
0. 30233 
0, 25000 
0. 11628 



0. 16071 
0. 07143 
0.50000 
0. 19643 
0. 07143 
0.0 



3 . 
0.2727? 
0. 18182 
0.30682 
0.14773 
0.02273 
0.06818 



23 
4 

0. 51020 
0, 04032 
0.32653 
0.0 
0.0 

0 122^5 



COLUMN TOTALS 



1.00000 



1* 00000 



1,00000 



! . 00000 



INTERDEPENDENCE TABLE 



I. AGRI-MINES 
2«C0NST-MANUF 

3. UTIL-C0MM 

4. HOUSEHOLDS 



2.00425 
0. 49906 
1. 32148 
0. 89918 



1 .49851 
1. 69^^46 
2 , 14468 
1, 1Q290 



1.51822 
0. 76318 
2, 8884? 
1.0356! 



1 5794R 
0. 5730? 
Ic 7049? 
1 , 



COLUMN TOTALS 



4.72396 



6. 441 55 



6.2 0 544 



5.69937 



INCOME MULTIPLIER 
TYPE I 

1 1.6147 

2 3 0483 

3 3,8059 



INCOME MULTIPLIER 
TYPE I I 

1 2,9742 

2 5.614H 

3 7 0103 

4 1.3419 



EMPLOYMENT MULTIPLIERS 



DIRECT EFFECT 

1 0*15698 

2 0.28571 

3 0 07386 



TOTAL EFFECT 
0 29238 
0. 55616 
0. 36747 



MULTIPLIER 
1*36256 
1.94658 
4. 97499 
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ENVIR. FACTORS PER DOLLAR OUTPUT 



2A 



ERIC 



1. SULFIDES 
2.800 

3. WATER REQ 



•0. 581395 
•0. 069767 
0. 767AA2 



0. 0 
-0.071429 
0, 785714 



■ 0 , 1 <^ 1 ^ « 

0.11 ->,f,-^h 
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ENVIRONMENTAL INTERDEPENDENCE MATRIX 



l.SULFIOES 

2. BOD 

3. WATER REO 



■0, 805911 
■0. 148953 
1.174313 



-0. 525704 
-0. 243712 
1. 6144 74 



■0^717265 
■0. 261B73 
1.049428 



ENV, SELF MULTIPLIER I 



I. SULFIDES 

2. BOD 

3. WATER REQ 



1.386167** ********** 



2. 134987 
1.530165 



3. 411963 
2. 054785 



3,Q44956 
2.711154 

9 2349 70 



ENVIRONMENTAL INTERDEPENDENCE MATRIX II 



l.SULFIOES 

2. BOD 

3. WATER REO 



■1, 478<330 
'0. 303120 
2 814453 



-1, 351''03 
-0.. 432308 
3 626207 



■0.4^'34^.1 
2 ,9 ^"^0^ 



-1 , ^ 7u..r, . 
1 -> a y 



ENV. SELF MULTIPLIER II 
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1. SULFIDES 

2. BOD . 

3. WATER REO 



1 2 
2.543758** *«*«**««♦» 
4,344728 6. 059314 
3.667317 4. 615172 



3 

8, ?0H: fr' ' ^ 3 ,1 ;•. (, 



7 ^ -J 



ENVIRON - EMPLOYMENT MULTIPLIERS 



1. SULFIDES 

2. BOO 

3. WATER REO 



■2. 756400 
■0. 509452 
4.016418 



-0.945231 
-0-438201 
2, 902871 



•1.951996 
•0.712636 

2:.855S1 5 



TYPE I ENVIRON - INCO^IE MULTIPLIERS 



I^SULFIDES 

2. BOD 

3. WATER REQ 



■1. 650880 
•0.305124 
2.405537 



-0, 877970 
-0,407019 
2 696308 



-1.275733 
-0.465770 

1,866523 



TYPE II ENVIRON - INCOME MULTIPIERS 



1. SULFIDES 

2. BOD 

3. WATER REQ 



! 

■1. 644750 
■0. 3 37107 
3,130014 



-1. 225133 
-0.392426 
3. 2S7876 



3 

•1,441 OPO 
•0.424^?? 

?« «37^7Q 



- 0 , 7 v P /. "7 c 



PROGRAMMER DOCUMENTATION 

The Input-Output model used in this program presents the economy of a 
particular region as a system of simultaneous linear equations. These equa- 
tions represent the distribution of each economic sector's output among its 
purchasers. Purchasers are separated into two main groups, the processing 
sector and final demand • The system of equations is solved to give the output 
of industries of the processing sector expressed as a function, of the final 
demand sector • The solution is in the form of coefficients that express the 
total change in output in each industry for a change in final demand in each 
industry. 

The environmental matrix is formed by expressing resource use in terms 
of units per dollar of output. When this matrix is multiplied times the 
above matrix of coefficients, the total change in environmental factors (posi- 
tive for resource inputs, negative for pollutants) per dollar change in final 
demand are estimated. Trade-offs between environmental factors and economic 
variables such as income and employment are estimated by combining relations 
found in the economic analysis with those derived in the environmental analysis. 

Basic Mathematical Model 
The following is a mathematical presentation of the model. 

Let, = x_^_. + ... x_^^ " 1,2,.. .n (1) 

where = total output of industry i 

x^j» purchase of output from industry i by industry j 

y « purchase of output from Industry i by final demand sectors. 
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It is assumed that a linear relationship exists between purchases of a sec- 
tor from other sectors and the level of output by that sf^ctor. Since total 
output equals total purchases, this may be oxprcssed as: 

X. . = a. .X. (2) 

or 

= ^J^'j (3) 

where, = total purchaiies by sector j. 

Substitutin-g (2) into (1) yields the ec'iation, 

X. - a. .X, + a_X^ c. X + y. ' 

1 ij 1 i2 2 ir. n -^1 (A) 

Data being fed into the progr:^'\ consists of purchases, x^^ of equation 

(1) , and is printed out in the data and r rv-isactions !^uK^'ix. The a., eleirn-ntb^ 

of equation (4) are calculated and present?.! as elements of vhe D!r».M::t Re- 

qiiirements r.^atrix, sometiraes referred to m Mi*: i.ochnUa.l rn^jf i c 1 triC::; of 

production. 

By re.irranging (A) (. his system ci. equc^ r.ion niay lie Xv'ritlren, 
1 i 1 12 13 J lu n • J 

a >: - (J - a,,,J:.., - a.^.^A -a., \ ^ v., 

. a -X, - a ^X^ - a ^X,;...-M. - a )X. - y 

nl 1 n2 2 n3 ' lui n n 

Equation (5) can then be expros:-,^.*! in matrj : notation as follows: 

(I - A)X = Y ) 

Inverting the matrix (I - A) yields 

X = (I - A)"'*^Y t7) 

The (I - A) ^ represents the matriA oi parcLril derivatives of X with respect 
to Y (i.e., ^X^/Oy^. )- Elenieuts of this nuitrix are commonly referred to as the 
matrix of in terdepeadence coefficients and represent the change in output of 
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industry i for a change In final demand in industry j . 

The environmental data is transformed into resource use (pollution) per 
unit of output. The relationship of resource use (pollution) and output is 
assumed to be a linear function also, 



r, . = b, .X, 



(8) 



where 

r, . - amount resource requirement or pollution output 
X_. = output of sector j 
and b- , = a constant. 

The technical coefficients b, ^ are thus derived from the combination of 
calculated output and collected resource data as follows, 

b, , = r- ./X, 
It j ^3 



(9) 



Let =^ total regional level (positive or negative) of environmental 
factor K. Thus, 



h ^hlh'-hlh^"-^ ^n^n 



^2 = ■^•••^ ^2n\ 



(10) 



R = b ,X, + b -X. +. . .+ b X 

m ml 1 m2 2 mn n 



Again, in matrix notation 
R = BX 



where R = 



m 



B = 



^1 ^2---^ln 



'21 ^22- ••^2n 



ml m2 mn 



X = 



X, 



X 



n 
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Substituting equation (7) from the input-output model into equation (!]) 
yields the environmental-output equation, 

R = B(I - A)'^^Y. (12) 
For simplicity, let P = B(I - A) "'"Y so that equation (12) becomes 

R - P(Y). (13) 

The P matrix contains the environmental-output multipliers which express 
the relation between a change in final demand in a given sector and the total 
change in a given environmental factor. For resource use (e.g., water, minerals, 
land) these multipliers are positive and may be interpreted as the total change 
in the use of the resource within tl^e region per dollar of final demand in a 
particular sector. For pollutants the multipliers are negative, but their 
interpretation is the same, A P^^^ element indicates the total change in 
environmental factor k per dollar change in final demand in industry j. 

Calculation of Multipliers 

Environmental-Income Multipliers 

The trade-off between environmental factors and income is found by using 
the total income effect and the resource-output multiplier. The total income 
effect is usually found as the first step in generating income multipliers. 
If households are exogenous to the processing sector then the total income 
effect is' 

(9) 2 = W(I - A)"^ 

where Z = a vector of multipliers indicating the change in wages in a 
sector for a dollar change in output for that sector (9W/3y), 
W *= the vector of wages paid per dollar of output. 

The type I environmental- Income multipliers are then found by dividing each 
environmental-output multiplier by the total income effect of the appropriate 
ERJC sector. 
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where Z. = the ith element of vector Z. 

J 

= matrix of environmental-incorae multipliers 

Wl\en the households sector is made endogenous to the processing sector^ 
the procedure is similar to the above. However, the vector Z is calculated 
directly in the inversion process that yields the interdependence coefficient 
By using the household row of the inverted matrix, we obtain the type II re- 
source income multipliers* The difference between type I multiplier as des- 
cribed above and type II multiplier is that the induced spending of household 
is included in the latter, 

Envirbnme n tal-Employmen t Mult ipliers 

Multipliers for environmental factors and employment interrelations are 
calculated in the same manner as the type I environmental-income multipliers. 
In place of the income vector Z, a vector of employment coefficients is used. 
These coefficients are found by dividing the total employment of a sector by 
that sector's total output. The resulting vector will then contain elements 
which are the average number of employees per dollar of output. Total change 
in employment for the economy is then estimated by the equation, 

(10) E = L(I - A)*^^ 

where E = the total effect on employment for a dollar change in each 
industry's output 

L = the vector of direct employment coefficients 
Then, environmental-employment interrelations are calculated from, 
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(11) R,, 




/E 



where 



a matrix of multipliers indicating the cliange in environ- 
mental factor usage for a unit change in each industry's 
employment. 



Environmental-Self Multipliers 

The relationship between the amount of environmental factors an industry 
uses ( or produces in the case of pollutants) and that which it causes to be 
used through its purchases of inputs may be similarily estimated from the 
model. This "environmental-self" multiplier is calculated by dividing the 
environmental output multiplier for a particular economic sector and envi- 
ronmental category by the corresponding sector's environmental factor use 
per dollar of output. 
That is 



(12) ^ 



/b, 



where 



a matrix of environmental multipliers relating the total 
change in regional use of factor k for each unit change in 
its use in sector j. 



PROGRAM DESCRIPTION 

The program is written in FORTRAN IV with single precision. The JCL 
cards presented in the user documentation are for a one step execution avail- 
able on tlie Texas A&M University data processing system. Other systems may 
require a two or three step execution and corresponding JCL, 

Routines 

MAIN: This routine reads data in, makes the calculations 
demonstrated in the mathematical, model, and prints 
out the results. 
INVERT: This subroutine uses Gaussian Elimination to invert 
the (I-A) matrix. The (I-A) matrix is placed in 
storage as DATA and passed to IKVl'^.RT along wiLli sup- 
port 1 n V ; I r tab 1 e s and v c t or s . 



In for mat i n n passed 

DATA 
L 
lER 
LL 
MM 



R eceived as 

X 
N 

lERR 
L 
M 



COMPRESS: Aggregation of collected data is done by this subroutine 
as an option. The main routine then does its calcula- 
tions with the reduced matrices and vectors. Four 
variables and one vector are passed to COMPRESS. Two 
matrices and one vector are placed in common with MAIN 
and COMPRESS and are accessed through BLOCK DATA. 
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Information Passed Received As 

NC NC 

NN NEW 

NO NO 

MN MN 

NVEC NV 

BLOCK DATA: holds matrices DATA and POL and vector EMP in coumion 

with MAIN and COMPRESS. DATA contains the (I-A)"*^ matrix, 
POL contains the environmental factors, and EMP is the 
vector of employment totals. 

Dimensions and Initial Data Statements 

MAIN 

The maximum size of the needed matrices and vectors are specified in the 
dimension statements, and selected ones are initially filled with zeroes by 
data statements. Matrices DATA and POL and vector EMP are placed on a common 
card for access by COMPRESS. 

INVERT 

Matrix X and vectors L and M are dimensioned initially for use in this 
subroutine. 

COMPRESS 

DATA is dimensioned as DATIX. The maximum size of matrix POL and vectors 
NV and EMP are given. A common statement connects DATIX, EMP, and POL with 
MAIN. 



BLOCK DATA 

The arrays held in common, DATA, POL, and EMP, are dimensioned and filled 
with zeroes. 



Should it be desired to increase or decrease the maximuin size of any 
matrix or vector, care should be taken to see that all dimension cards 
containing the array have been properly altered. Also that, if applicable 
the proper number of zeroes are placed in the data statement. 



ERIC 



FJ.aW CHART 



START 







CONTROL 


CARDS 




f 


'^^AD 




ECON CARD 

' , — d 




CALCULATE 

t»RODUCTION 
COEFFICIENTS 



CALL 
INVERT 



1 



MAKE H 
INDOGEN 


OUSEHOLD 
OUS 








CALCULATE INCOME 
MULTIPLIERS 




NO 



CALCUIAT t 

EMPLOYMENT 

MULTIPLIERS 




CALCULATE 
ENVIRONMENTAL 
COEFFICIE NTS 



CALCULATE 
ENVIRONMENTAL 
INTERDEPENDENCE 
MATRIX (TYPE 1) 



CALCULATE TYPE 1 
ENVIRONMENTAL "SELF " 
MULTILIERS 



1 

CALCULATE TYPE 11 
ENVIRONMENTAL 
INTERDEPENDENCE 
MATRIX 



V 



CALCULATE TYPE 11 
ENVIRONMENTAL "SELF 
MULTIPLIERS 



CALCULATE 
ENVIRONMENTAL - LABOR 
MULTIPLIERS 




NO 



1 



CALCULATE TYPE 


1 & 1 1 


ENVIRONMENTAL - 


INCOME 


MULTIPLIERS 






PROGRAM LISTING 
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FORTRAN IV G LEVEL 21 



MAI 'J 



OATF 



7-> 1 7 



1 0 / ' ' 



0001 
0002 



0003 
0004 
.0005 



0006 
0007 
0008 
0009 
0010 
0011 
J012 
0013 
0014 
0015 
0016 
001 7 
0018 



0019 
002 0 
0021 
0022 
002 3 
0024 
0025 
002 6 
002 7 
002 8 



0029 
0030 
0031 



0032 
0033 
^ 0034 
ERIC 5 



c 
c 
c 



c . 
c 
c 
c 

c 

10 
15 



20 

c 
c 
c 
c 
c 

25 

?0 
C 

c 
c 
c 
c 

c 

35 
40 
45 



REAL MUL1,MUL2 

DIMENSION OATA( 100, lOOJ f EMP{ 100) , POL (2 5, TOO) , 
1NAME( 100,4 > ,NAM2{25f4> ,CTOT ( 1 00 ) , RT^T ( I 00 ) , 
2DINT(90,90) ,PINT(25 ,90) ,DTEC(90,90> ,STOt^ (?5,'5S) , 
3MM( 100) , NVECJ 95 ) , NAME 2 ( 1 , 4 I , S TOR ? ( 100) ,M!JLl ( 1 00 > , 
411(4), JJ(4»,LL{1 00) ,DAT2(95,95) , MUL ' ( 1 00 ) , X ( 4 ) 
COMMON /OEP/ DATA, EMP, POL 
DATA I I /4*0/, JJ/4*0/,X/4*l 



lNAME(100,4»,NAM2{25f4) ,CTO 
2DINT(90,90) ,PINT(25 ,90) ,DT 
3MM( 100) tNVECJ 95 ) ,NAMF2(1,4 

00 ) , X( 4) 

COMMON /UtP/ UATA,EMP,pnL 

DATA II /4*0/, JJ/4*0/,X/4*0.0/ ,NVt^C/9S*0/ 
DATA IER/0/ 

« « « « « « 4c « « « « « 4c « « 4c )(c « « « lOc 4t 

* CONTROL CARD RLAD * 

4c4t4c**4ci(<i(t**<<4e4c4i4i4e>(t*4c 

READ (5,5) NN,NP,NV,NI 
READ (5,5) MI>l,LAB,NCOMP 
FORMAT (615) 
MNA1=MN+1 
NVA1=NV+1 
NNAl =NN+1 
NNA2=NN+2 
NNA3=NN+3 
NNA4=NN+4 
NNA5=NN+5 

NPAl = NP'rl 

NPA2=NP+2 
NPA3=NP+3 

4c 4c 4c « 4i « « 4t « 4i * « 4c 4t « « « * 4i 4t 4c 4c « i«c 4e ft « 

* DATA MATRIX IS RFAC IN * 

4c * 4c 4c * 4< 4c 4c 4c « 4c 4c « « 4c 4< 41 4t 4< 41 41 4c 4< 4c 4c 41 « 4c 

READ (5,15) { II ( 1) , JJ ( I) ,X( I) , 1 = 1 ,4 ) 
F0RMAT(4(2I3,F12.3) ) 
IF (11(1. ).E0.-99) GO TO 25 
00 20 J=l,4 
K=I I ( J) 
L=JJ( J) 

IF ( K.EQ.O.OR.L.EO. 0) GO TO 20 
DATA (K,L)=X(J) 
CONT INUE 
GO TO 10 

4e 4c 4c V 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c « 4c 4i 4c 4< 4c 4c 4c 4c 4c 4c 

* EMPLCYMEr^iT VECTOR READ * 

4c 4c 4c 4c lie « 4c 4c 4c 4c « 4c 4c 4c 4c ^ 4c 4c 4c 4i 4c 4c 4c 4c 4c 4c 4c 

IF (LAB.EQ.O) GO TO 35 
READ (5,30) (EMP<I ) ,I=1,LAB) 
FORMAT (6F10.0) 



4e4c4>4c4c4c4c4<4>4c4i4c4c4c4c4c4c4i4c4i4c4c4c4>4c4<4<4c« 4c4c 

* RES MATRIX IS READ IN * 

vc 4c 4c 4c 4< 4c 4c 4i 4c 4c 4c 4c 4c 4c 4< 4c 4c 4c 4c 4< 4c * 4> 4c 4c 4c 4c 4i 4< 4c 4c 

IF {MN. EO.Ol GO TO 55 

READ (5,45) I 1 1 ( I) , J J ( I ) , X U ) , T = 1 ,4 ) 

FORMAT (4(213, F12. 2) ) 

IF ( I I ( 1UEQ.-99) GO TO 55 
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MAIN 



0 A T F = 7 117"' 



0036 
0037 
0038 
0039 
OOAO 
OOAl 
0042 



0043 
0044 
0045 
0046 
004 7 
304 8 
0049 



0050 
0051 
0052 
0053 
0054 
0055 



0056 
0057 

005 8 
0059 

006 0 
0061 
0062 
0063 



ERIC 



0064 
0065 
0066 
0067 
0068 
0069 
0070 
0071 



50 
C 

c 
c 
c 

c 

55 
60 
65 



70 
75 
C 
C 

c 
c 
c 



RO 

85 

C 

c 
c 
c 
c 



90 

95 
100 

C 
C 

c 
c 
c 
c 
c 



10 5 
110 



DO 50 J=l,4 

IF ( TI( J».E0.0.0R.JJ(J).EQ.0> GO to SO 

K=II ( J) 
L=JJ { J) 
POL{ K,L )=X{ J) 

CONT INUE 
GO 40 

* l^EAD REGIONS AND TITLES * 
DO 6 0 I=1,NN 

READ (5,65> ( NA ME ( I , J ) , J =1 , 4? 

FORMAT {I3f3X*3A4) 

IF { m, EQ. 0) on TO 75 

00 7C 1=1, MN 

READ (5,65) ( Nfi M2 { I , J ) , J=l ,4 ) 
CONT INUF 

« CALL SUBROUTINE COMPR'^SS * 

IP { NCOVP. FQ. 0) GO TO 85 

NC=NC0MP*3 

READ (5,80) NO, (NVFCd) ,I = l,NCI 

FORMAT (1315/, (1215)) 

CALL CMPRSS (NC,NN,NO,MN,NVFC ) 

CONTINUE 

:tc 4ci«c 4c ft « i«c i«c « i«c 4< itt « itc « 

* DATA MATRIX IS PRINTED * 

DO 95 I =1,NN,9 
K = I + 8 

IF (K.GE.NN) K=NN 
WRITE (6,90) (L,L=I,K) 
FORMAT ( '1' ,14X,9(10X,I2) > 
DO 95 J=1,NN 

WRITE (6,100) ( NAME( J ,L ) ,L = 1, 4) , (DATA { J ,.M) ,N= I ,K ) 
FORMAT ( I4,«, • ,3A4,2X,9F12.3) 

« « 4c ^ * :(c « << « « :4c :(c « « 3(t * 4c « 4c « « * 4c ♦ « 4t * * a >Cc « « >Cc :«( 4c >^ * « « * # * « « * >!: 4t « >k 'tc >!c 

4c TOTALS G SUBTOTALS OF COLUMNS C ROWS * 

* CALCULATE PRINT & * 
« STORE * 

4c4c4c4c4i4c4c«ft4c4c4c4c4i4c4c4cA*4[4c4c4c4c4c4c4c4c4'4c«4c4c4c*4c4c4c4c4c«a!c*i!c«*4ri(r«>v^, 4i4r 
SUMR=0. 0 

SUMC=0.0 

DO 110 1=1, NN 

RTOTd )=0.0 

00 105 J=l,NN 

RTQT ( I ) =RTrT{ I ) ♦DATA( I , J) 

SUMR=SUMR+RTOT( 1) 

00 120 1=1, NN 
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MA IN 



DA- 



0072 
007 3 
0074 
0075 
0076 
Q077 
0078 
0079 
0080 
0081 
00fl2 
0083 
0084 
0085 
0086 
0087 



008 8 
0089 

009 0 
0091 
0092 
0093 
0094 
0095 
0096 
0097 
009 6 
0099 
0100 
0101 
0102 
0103 
0104 
0105 
0106 
0107 
0108 
0109 
Oil 0 
0111 
0112 



0113 
0114 
0115 
0116 
O 117 

ERIC 11 8 



115 
120 



121 
122 

123 
124 

12 5 
C 

c 
c 
c 
c 



12 c 
130 

13 5 
140 
145 
150 

155 
160 

165 



170 

C 

C 

C 

C 

C 

c 



CTOT( n=o.o 

DO 115 J=1,NN 

CTnT{ I)=CTOT( I >«-OATA( J,I » 

SU^^C =SUI«C + CTOT{ r » 

DO 121 1=1, NN 

DO 121 J=1,NP 

DATA (I ,-NNAl) =DATA( I ♦ NN Al ) +0A T A (I , J ) 

DO 122 1=1, NN 

DATA ( I ,NNA2» =RTOT( n 

00 125 I=1,NNA2 

DO 123 J=1,NP 

DATA (I ,NNA3) =nAT A( I , NNA? ) +DATA( J , ! > 
00 124 K = .'vJPA] ,NV 

DATA ( I ,NNA4) =DATA( I , NNA4 J + DA ta ( K , I ) 
DO 125 !.=NVA1 

DATA ( r,NNA5) =0ATA( I , NNA5 ) + DA T A < L , I > 

* TRANSACTIONS MATRIX IS PRINTED * 

00 160 I=1,NNA2,9 
K=I + 8 

IF (K.GT.NNA2) K=NNA2 
WRITE <6,126) {L,L=I,K) 
DO 150 J=1,NN 

WRITE (6,145J ( NAME( J ,L) ,L = 1, 'f) , (DATA( J , N ) = 



I , K1 



IF ( J.rO.NP) 
IF (J.Cg.NV) 

N1=NP 

IF ( J.FO.NI ) 
IF (J.fU.NI) 



WRITE (6,130) (DATA(L ,NNA-^ 1 ,L = 1 ,K ) 
WRITE (6,135) <nATA{L,NNA4) ,1 ^ ( ,K > 



(DATA(L fNNA'i ) ,1^ T ,^ i 



SUM =',f]7,C] 



175 



W^ITF (6,140) 
GO TO 150 

FORMAT! '1 TRANACTION TABLE'///' • , \ X , « ( 1 0 X , ' ^ M 
format; //2X,' INT, PURCHASFS • , 3 X , 9F 1 2 . 2 / / ) 
FORMAT ( //2X, • VALUE CREATED • , 3 X , 9 F 1 2 . ?/ / I 
FORMAT( //2X TOTAL IMPORTS • , 3 X , 9 F 1 2 . 2/ / » 
FORMAT (14, ' ,3A4,2X,9F12.2» 
CONTINUF 

WRITE (6,155) (CT0T(M4) ,M4 = I, K) 

FORMAT( //2X , 'COLUMN TOTALS ' ,3 X, 9F1 2. 2// I 

CONT INUE 

WR ITE (6,165) SUMR, SUMC 

FORMAT! '1',' ROW SUM = • , F 17 , 4 , 5X , • C OL U'^N 
DO 170 1=1, NN 
DATA (NNA1,I )=CTOT(n 
CONT INUF 

* ENVIRONMENTAL MA"^RIX G EMPLOYMENT * 

* VECTOR IS PRINTED * 

4i itc « 4r 4i 4i * i4i ««« 4ii{c « i» iti Id :4i ** 4c 4i i4c]4i 4: i«c 4c 



IF ( MN. EQ.O) GO TC 185 
DO 180 I=l,NPAl,9 
K=I + 8 

IF (K.GT.NPAl) K=NPA1 
WRITE (6,175) {L,L=I,K) 
FORMAT! '1' ,3X,' ENVIRONMFNTfiL 



FACTORS ', /////I 5X, 



FORTRAN IV G LEVEL 21 



MA IN' 



0119 
0120 
0121 
0122 
0123 
0124 
0125 

012 6 



0127 
012 8 
012 9 



0130 
0131 
0132 
0133 
0134 
0135 
0136 
013 7 
0138 
0139 
0140 
0141 
0142 
0143 
0144 
0145 
0146 
0147 
0148 
0149 
0150 
0151 



0152 
0153 
0154 
0155 
0156 
0157 
0158 
0159 
0160 

ERlO^i 



IPO 
1B5 



190 

195 

C 

C 



200 

C 

C 



c 
c 

205 



210 



21 5 

220 
225 
230 
235 



240 
245 
C 
C 

c 
c 
c 



246 



19(10X,I2)) 
DO 180 J=1,MN 

WRITE (6.145) ( NAV2 ( J ,L » ,L =1, 4» , (P£)L( J,N t ,V'= T , K » 

CONTINUE 

CONTINUE 

IF (LAB.EQ.O) GO TO 195 

WRITE {6fl90) { I,exP(n,T = l,t^P| 

FORMAT { n ' f 5X t • EMPLOYVENT VE CT OR • , / / , ( S X , I 7 , • , • , - v , 
1F10.2 » ) 

CONT INUF 

DATA MATRIX IS STORED AS 0AT2 

00 200 I=1.NNA1 
00 200 J=1,NNA1 
DAT2( I , J) =OATA{ I , J) 

* DIRECT REQUIREMENTS ( I £2) * 

* CALCULATE e PRINT * 

CONTINUE 
00 210 !=1,N1 
CTOTC I )=0.0 
DO 210 J=l f NN 

DATA (J, I )=DATA( J,I ) /D ATA ( NNAl , ! ) 
CTOT( I ) =CTOTt n +DATA( J ,n 
DO 230 I =1 f NP Al ,Q 
K=I+fi . 

IF (K.GT»N1» K = N1 
M2 = K 

WRITE (6t215> (L,L=ItK) 

FORMATCl DIRECT W E Q!) I R E M FNT S ' / 1 '^X , Q (i ) x , T ? M 
DO 220 J=1,NN 

WRITE <f5,22 5) (NAME(J,L),L = 1,4) ,(nATA(J,N) ,N = T ,Kt 

FORMAT (14, « ,3A4,2X,9F12«5) 

WRITE (6»235> ( CTCT { M4 ) , V4 = I , k ) 

FORMATC //2X, 'COLUMN TOT ALS • , ''i X , 9F 1 2 . » , 

IF (Nl.EQ.NP) GO TO 245 

00 240 1 = 1. Nl 

DO 240 J=1.N1 

OTeC( I, J)=OATA( I , J) 

CONT INUE 

* INVERT MATRICES 1 £ 2 * 

L=N1 
M=l 

DO 246 1=1. Nl 

DATA ( I .1 )=1,0-DATA( I , n 

DO 246 J=1.N1 

IFd.NE.JI DATA(I,J»= -OATA(I,J) 
CALL INVERT ( DATA , L . I ER t LL , MM » 
IF ( lEP.NB.O) STOP 
DO 2 50 I =1 .Nl 
CTOTt I » =0.0 
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MAIN 



^177 



0162 




DO 250 J=l ,N1 


0163 


250 
C 


CT3T( I ) =CTOT( I ) +DATA{ J, I ) 




c 






c 


* PRINT L STORE INVERSES 1 2 * 




c 
c 




0164 




IF (Nl.Eg^NPAl) GO TO 260 


01 65 




DO 2 ^5 J =1 f Nl 


0166 




DO 255 J=lfNl 


0167 


25 5 


0 INT( I , J » = r3ATA( I , 


016 8 


260 


CONTINUE 


0169 




DO 275 I=ltNl »9 


01 70 




K = I^8 


0171 




IF (K.GT.Nl) K=N1 


01 72 




WR ITE { 6,265) ( Lt L=I , K) 


01 73 


265 


PORMAT( n INTERDEPFNDENCE T A BL F • / 1 5 X , 9 P OX , ! ) ) 


01 74 




DO 270 J=l ,N1 


01 75 


270 


WRITE t 6,22 5) ( NA ( J , L ) t L = 1 , 4 ) , (OATA( J,M^ ) , m4= t,k ) 


0176 




WRITE ( 6,235) ( C TOT ( M4 ) , ^^4 = 1 , K ) 


0177 


275 


CONTINUF 


0178 




N1 = N1 + 1 


0179 




IF (N1»EQ.NPA2) GO TO 28b- 


01 80 




00 276 I=1,NNA1 


01 81 




DO 276 J=1,NNA1 


0182 


276 


DATA ( I ,J)=0AT2( I , Jl 


01 83 


C 


GO TO 205 




c 






c 


'0^ AT THIS POINT THE STORED MATRICFS ARF ^S FH^^l/^W^ ^ 




c 


* DINT=INVERSEI nAT2=nRIGrcnNnATA * 




c 


* DTEC=TECH COEFF W/ HH POL =ORIG ^^FS OA^^ * 




c 


'0^ OAT A = INVERSE I I * 




c 






c 
c 






c 


* I NCOME MULT I PLIERS I &I I « 




c 




01 84 


c 

28 0 


DO 290 1 =1 ,NP 


0185 




VALUE=0«0 


0186 




DO 2 85 K=l , NP 


0187 


2 8 5 


VALUE=VALUE+DINT(K, I )*0TEC(NPA1 ,K) 


01 88 


290 


ST0R2( I ) =VALUE 




C 


TYPF I 




C 


CALCULATE & STnoc 


0189 




DO 295 1=1, NP 


01 90 


29 5 


MULl < I ) =( ST0R2( I ) /DTEC(NPA 1,1 ) ) 




C 


TYPF II 




C 


CALCULATE L STnoc 


0191 




00 300 I=1,NPA1 


0192 


300 


MUL2 (I )=(DATA(NPA1 , I)/DTFC(NPA1, I ) ) 




C 






C 




0193 




WRITE (6,305) ( I , MULl (I ) , I =1, NP ) 


0194 


305 


FORMAT! ,lbX,» INCOME MULTI PL I FP • / • •,^^X,»'^Yor 






l/(6X,I2t6XfF9*4)) 


0195 




WRITE (6,310) (I ,^UL2 ( I) , 1 = 1, NPAl J 



01 96 



019 7 
01.9 8 
01Q9 
3? 00 
0?01 
020? , 

020 3 
0204 

02 06 
0207 

0208 
0209 



0210 
0211 
0212 



310 
r 

C 

c 

G 
r 



n 5 



?2 C 

^^0 
3^5 



336 
34C 
r 



r 

r 



Fn = MAT( » 1 * ,1 OX t M^COMF MULTTPL I /• '^isy.'Tvor rt', 
1 /( 6X, 12 , 6X,P9.4) ) 

A FMPLOYMPNiT MIJLTTPLIFPS * 
* C^LCULATr:, P" I NT , f, STOP= IN P AT a ( , n *■ " 

IF { LARo FO. 0» GO TP "J^^O 
DO 315 I=I,NP 

F Mp( n = Fwp( J > /OAT? ( NVAl, T ) 

nn 3 25 T =1. ,NjD 

CTOT ( I ) =0o 0 

VAL'JF = D.O 

nn 320 K=] ,NP 

VAL!JF = VALUE+nriT(K,n*CMP(K) ' 

TT IT { I ) =VALUr 

no 3 30 J=1 ,NP 

RTIT { J) =CTnT( J) /fc-MP( J ) 

WPITP ( 6,335» { I .i^'-^PC T ),CTnT( I ) ,cT"i^ ( T ) , ! = ■' 

FC'?MAT{ • I' ,2BX,•F^'^L'^YMFNT wiJL T [ PL T 'P '">'//' ^ X, 
I'DI^FCT f^FFcCT' ,9X, 'TOTAL ^FF t^CT • , OK , • II. I PL t r c i , 
2/( 10X,I 2,3X,F c, 5,1 -^X.Fn, 5,10X,Fo„c;n 

nn 336 T=l ,NP 

DATA {NNA4 , 1) ^RTP"^ ( T ) 

COMTINiJc - - ■ • 

itc ;tc 4c « 1^ itc :«c >!t # « >!e « 3}c i{: « )te 1^ >!t # :«c itc !«( itc 4c « « # )te « « i)c :4e :tc :(c >!t ift 

*- ^NVT?; ONMENT AL FArrrips pfo i nuTDijr * 
*CiVLCUL4TE PRIMT F. ^TORE IN PHLKItl) 

14( }tc « )tc « « 4: )tc )tc 4c >!c <(:tc « )Sc * <( « 4: 4c t Ik « iC 3(r » 4: « lit <!t ^ * * 



021 3 
021 A 
0215 
021 6 

021 7 
021 
0219 
0220 
0??l 
0222 
0223 
0224 

022 5 
0226 
0227 



I 022 8 
f 0229 
i O 12 30 

;ERIC23i 



34 5 



350 
35 5 
360 

365 

C 
r 

c- 
c 
c 
c 



IF ( m, EQ.O) r,r TP s?^ 
DO 345 T =1 ,MPA1 
nn 345 J=1 »MN 

POU JtT )=PnL{.J,I > /DAT 7 < MM A 1,1 ) , 
DO 3 65 I-1,^JPA^ ,9 
K = I+ S 

IF (KiGT.NPAll K = NPai 

IF (I/2*2,FQ.II WRITF (6t350> (l,L=>,K) 
IF {I/2*2,cQ, T» GO TP 360 
FORMAT (/////15X,9(10X,I?) t 
WOITc (6,355) (L,L=I,K) 

FQPM AT( • IFNIVI FACTORS PER HOLLAR lyTr*' J ^ ' / r / 1 5 X , 
l9(10X,I2n 
CONTINUF 
no 365 J=l ,MN 

WRITF (6,390) < NA"? ( J , L ) ♦ L =1 , 4 ) , ( P^L ( J , » , ' 1= T , K » 

* 4c 3tc 4c 4c 4i 4c 4i 4c * 4c 4c 4< 4c 4< 4> 4< 4< * 4i * 4 4< « 4c 4i 4c # 4c 4c A 4r ft ic 4c 4c 4c A itc '«c A >!c 

* ENVIRONMENTAL T MTFR DFPFNOFNCF MATPIX - 

* CALCULATP, PCTNT,C STQCF IM oimt 

4c 4c 4c « 4: 4c 4c 4c :^ 4c 4c >tc 4c 4c 4c 4c 4c 4i « 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4c 4: ft A.4c 4c 4c 4c 4: 4e 4c « 

DO .3 70 I =1 ,MN . 

DO ^370 J = l ,NP ) 

prn{i,j}=o 

on 3 70 K = l ,NP 
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TAT - 7 « 



2 0/2 



023? 
02^3 
0234 
0?3.'5 
02 3 6 
-.02>7 
023R 
0?^9 

0240 
0241 
0242 
0243 
0244 



^70 



37'> 

3fl5 

395 
r 

C 
C 
r 

C 

r . 

r 



P TNT ( I , J ) =P INT ( I ♦ J 1 + POL ( I , K ) * DT ^■ 
nr* 39^ T=lx,NP,<? 
K = T+« 

IF (K.GT.NPI K=NP 
IF v( !/2*2.FQ. I ) WRTT=. {6,350> 
IF ( I/2*2,FQ. T I r,r TO 3B0 



{K , J » 



{L ,1 



WPITE (6,375) (L,L^ 



FG^MAK M i^MVinnMvrtvjTAL T NTPr nfPFNriPrgr ' 
1 'MATPIX'///15X,Q{10X,I?) ) 
CnNTlNt.iF 
no 3f^5 J = 1,MN 

WPTTF (6,390) ( NAW? ( J ,L) ,L=1 
F (If MAT (lA, 'o • ,3A4,?X,9Fl2,f ) 

rnr^JT TMUF 

* FNVIR, wtiLTTqiTFOS * 

* CALClILATc PC INT * 



^ Y P"^ T 



0245 
0246 
024 7 
0248 
0249 

0250 
0251 
0252 
0253 
0254 
0255 
0256 
02 57 
0258 
0259 



0260 
0261 
0262 
0263 
0264 
0265 
0266 
0>67 
02 6 8 
02 69 
0270 
0271 



40 0 



40 5 
410 

41 5 
C 



42 0 



42 5 
430 



nr 400 i=i,MN 
nn 400 J=l ,NP 

IF ( POL ( I , J)«Nc,0< 



0) nAT2 ( I , J ) = ( PT^iT ( T , J y/p-^i ( T , J n 



'^^TJ[ I , J I =0, 0 

HAT ■) ( ! , .1 1 



TF ( POL ( T,JK i^O. .:>. O.&MO. PINT ( [ , J I , FQ . 0 , ■) ) 
IF (POL( T ,J).FfJ3 0, O^.'NrJOo PINT(T,J).Nr.O^M 
1=10 30000 
rn 4 15 T=1,ND,9 
K = I+3 

IF (K.GT.NP) K=NP 

IF ( 1/2*2. FO. n W^lTr (6 ,350) (L,L = I,K) 

IF ( i/?*2,'^Q, T ) r,r TO 41 0 

WRTTF ( 6,405) (L,L = T,K» 

FQt'MAK'l ='gv. SPLF yui.TIPLIFP I • / / / Y , ( ITX , I 7 ) I 
COJ^TINUF 
OD 415" J = 1 ,MN 
WRI TE ( 6, 390) 



S0 



( NAWo ( J ,L» ,L=1 ,4) , (DAT?< J,H) 



73 



TYPF I T 

FNVIPflNMFNT I ^jt cr< h= P c\i ncMf p 
MATQIX STO'?- p.' oTpr 

on 4 20 .I = 1,MN 
00 420 J=1,NPA1 
DT=C(I,J)=0. 
nn 420 K=l ,NPA1 

OTEC ( I , J)=OTFC( 1, J) + oPL( T f K»*nATA(K, J) 

on 440 1=1 ,NP&1 ,9 . • 

K = \+^ . .- 

IF (K.GT.NPAl) K=NPA1 

IP (I/2*2.P0.I) WPTTF (6,350) (L,L=I,K) 
IF ( 1/2*2. FO. n r,G TO 430 
WRITE (6,425) (L,L=I ,K) 

FORMAK'l FNVIROMMFNJT AL I NTB ROFPFNOFNC - 
1 'MATRIX II •///15X,Q(10X,I2) ) 
CONTINUF 
DO 4 35 J=l , VN 



MAIM 



?o/; 



0?75 



0276 
0?77 
0778 
0?7Q 
0?flO 

02R1 
02 fl? 
02P3 

02fl5 
0? fl6 
0287 
02 8 8 
0289 
0290 



0291 
0292 
0293 
02 9« 
0295 
0296 
0297 
0298 
02 99 
0^00 
0301 

0302 
0303 
0304 
0305 



0306 
0307 
. 0308 
03 09 
0310 
0311 
0 312 
O 13 

ERICU 



43 5 
r 

c 
c 
r. 



44 5 



450 

45 5 

46 0 
r 

c 

c 

c ■ 

.c 

c 

c 

c 



46*^ 



470 
'4 7 5 

480 

48'^ 

C 

C 

490 



WPITP ( 6,390> ( \/lv?( J,L) tL»lt4) t (OTFr<j,.M) ,v = T,x'> 

▼ Y f r T T 



no 445 I=ltMN 

DO 445 j = l ,MPAl 

IF (DOLU , J». NF , 0.0) Pf^L < I ,J» = n"^<^r( T , J) / Pf 

TF ( POL ( I , J) ,FQ, 0,0o AMO, DTEC( T , J ) , =Q , 0 , 1 

IF ( PnL( I , JUFQ.vO.Oo AND. OT "^r. ( T , J ) , NF , 3 . 1 i 



L ( T , J > 



llOOTOOO 
pri 4 I =1 t.NnAl , " 
K"= T + 8 

IF ( K.GT.NPAl ) K = NPAV 

IF ( I/2*2<.rQ. T) WPI"^'^ (6t350) (1,1. = t,k) 

IF ( 1/2*2. fo, T » nn -^n 4'=;5 

Wa ITF ( 6,450) { L f L = I tK, ) 

FH^MAK'l Ff^jv, S^LF f^ULTIPLlFo I 1 ♦ / / / 1 X , ^ ( " 0>: , I ? > ). 

CONTINUE 

00 460 J=l ,MN 

WPITP { 6,390» ( MAW?( JtL) tL=l,4J , (PnL(J,'H ,N=!,'<) 

Xt *****«■(* !fr ************ * 

« CO VBI NATION MI.ILTTPL IFPS * 

* rfyjvi PONMCMT/VL _ FMPLOY * 

* FNVIRnNVFNTAL - TNCOMF * 
•*******>)<*******«*>(<**>■«*****♦**** 

FNVIC OMMFWT AL _ rMO| nyvrM-r 

r LC ','1. r t'. TMT 

IF ( LAB.EO.OJ r,n Tfi ^.fi5 
no 465 J=l ,NP 
Qf^ 465 Islt^'M 

STOR ( I , J) =PIN-^( I , J» /TTIK J ) 
on 4fl0 T=l ,NP,Q 
K = i + «( 

TF (K.GT.NP) K=NP 

IF ( I/2'!«2,F I) WOT Tc . (f,,350) (L,L=T,K» 
IF ( 1/2*2. C-Q.T) Gr TO 475- 
WRITE ( 6,470» (L,L = I,K> 

FnoMAT('l ^^NVI^^CN - c-MPLnyvifMT vi )LT I PL I^ "^^ ' / / / I X , 
"l9{10X,I2)) 



CO'VTtMUE 

00 480 J=l ,MN 

WP I TF ( 6,390) 

CGNTINUF 



( N^w' ( J»L) ,L=1,4) , ISTOf { J,M) ,M = T , K ) 



CNVIP ONMFNT At - T^^nvfr j 

C ALCUL-*"^ = 'I P"^ PIT 

on 490 j=i>NJp 

00 4-90 I =1 ."^N 

PIMTd t J» =PIMT( I ,J)/ST^R?< J) . 
on =i05 I =1 ,NP,9 . . 
K = I*8 

IF (K.GT.NP) K=NP 

IF (I/2*2.Fg.I) WPITF (6,350) a,L = I,K) 
IF .( 1/2*2. Fg, T) or TO SO.O 
WRITE (6,^95) ( L,L=I ,K» ■ 



0^1 7 



0319 
0320 
03?l 
0322 
0323 
0324 
032 5 
0326 
0327 
03?B 

0329 
0?30 
0331 
0332 
0333 
0334 



IV G LFVPL 
495 

50C 

r 
C 
C 



21 



MAIN 



? c/r 



19(10X,I?)> 
nn 5 05 J=l,yN 

WPIT^ (5,390) '(^JAM7(J,LJ»L = lt4), (OI^!T<J,^^,^;rT,K J 

FNVIT)M*^ = ''!'^ AL - tnithmc tj 

nn 510 J=1,NPA1 
no 510 1=1, MM 
510 aTFC( i , J»=OTFC( I , J) /nATA(NPAl ,J) 

no 525 i=i,NjDaT 

K = I*-8 

IF (K.GT.NPAl) K=NPAi 

IF (1/2*2. FQ.n WPTTTF (6,350) (L,L = T,KI 

IF ( 1/2*2. PQ. I ) r,C Tn 570 

WPITF (6,515) (L,L=T,K) 
515 FPRMAKM TYPF TI cvjvTRnN - INCOMF It T P ! 1= " S • / / / ""^ X , 

19(I0X,I2)) 
520 COMFINUF 

no 525 J=l ,MN 

WPITE (6,390) ( NAM7 ( J ,L) ,L=1 ,^) , (OTFC( J, ,M = '' t 
525 CQNTINUF 

STHP 

FMO 



t 
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INVERT 



0001 



0002 

0003 
0004 
00 05 
0006 
0007 
0008 
0009 
3010 
OOIl 
J012 
0013 

0014 
0015 
0016 
0017 
0018 
0019 

0020 
0021 
0022 
0023 
0024 
002 5 

0026 
002 7 
002 8 

002 9 

0030 
0031 
00,32 
0033 
0034 

003 5 

)36 



C 
G 
C 
C 

c 
c 
c 

c 
c 
c 
c 
c 
c 
c 
c 



SUBROUTINE I NVE RT ( X , N , I ERR ,L , M ) 

REAL MATRIX INVERSION SUBROUTINE PY GAUSSIAN 

ELIMINATION 

THE FOLLOWING CALLING SEQUENCE SHOULD BF USEH rn fNTfc 
THE FOLLOWING CALLING SEQUENCE SHHULO BE USEH TC ^MTpp 
THIS SUBROUTINE ..... X 

CALL INVERT<X,N', IERR> 

WHERE 

X IS THE MATRIX TO BE INVERTED 
N IS THE ORDER OF X 

lERR IS AN ERROR FLAG OENOTING SUCCESSFUL 
NON-SUCCESSFUL INVERSION OF X • 



DIMENSION XdOClOl) ,L(100),M(100) 
SEARCH FOR LARGEST ELEMENT IN X 
DO 80 K = 1,N 
L(K) = K 
M(K) = K 
XBIG = X(K,K) 
00 20 I = KfN 
DO 20 J = K,N 

IF(ABS(XBIG»,GE.ABS(X(I,J) n GO TO 20 
XBIG = X(( I , j) 
L(K) =1 
M(K) = J 
20 CONTINUE 
; NOW INTERCHANGE ROWS 

IROW = L(K) 

IF(L(K).LE.K) GO TO 35 

00 30 I = 1,N i 

WAIT = -X(K,n *~ 

X(Kf I) = X( IROW,n 
3 0 X(IROW,I ) = WAIT 

NOW INTERCHANGE COLUMNS 
35 ICOL = M(K) 

IFJ_M(K).LE.K) GO TO 45 

DO 40 J = 1,N 

WAIT = -X(J,K) 

X(J,K» = X(J,TCOL) 
40 X( J, ICOL) = WAIT 

NOW DIVIDE COLUMN BY MINUS PIVOT 
45 00 55 IC = 1,N 

IFdC.EQ.K) GO TO 55 

X(ICfK) = X( ICf K»/(-X(K,K) ) 
55 CONTINUE 
; NOW REDUCE MATRIX 

DO 65 I = 1,N 

DO 65 J = 1,N 

IF(I.EQ.K) GO' TO 65 
-IF( J.EQ.K) go" TO 65. 

X(I, J) = X(I,K)*X{K,J) ♦ X(I,J) 
65 CONTINUE 

: NOW DIVIDE ROW BY PIVOT 

DO 75 JRO = 1,N . 
IF( JRO.EQ.K) GO TO 75 
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INVERT 



I) A T P = 



0038 
003 9 



0040 
OOAl 
0042 
0043 
0044 
0045 
0O46 

0047 

0048 

0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 

0058 

0059 

006 0 

0061 

0062 

006 3 

0064 

006 5 

0066 



X(K,JRO) = X(K, JPC)/X(K,K) 

75 CONTINUE 

C CONTINUE PRODUCT OF PIVOTS AND PFPLACE oivOT f^Y 

C RECIPROCAL 

CALL 0VERFL(J1> 

CALL DVCHK{J2) 
.X(K,K) = I,0/X(K,K) 

CALL DVCHK(J2» 

IF( J2..EQ. n GO TO 600 

CALL OVERFL( Jl) 

IF(Jl.EQ.l» GO TO 600 
C NOW CONTTnUE OVERALL OPERATION 

80 CONTINUE 
C NOW FOR FINAL ROW AND COLUMN INTEi^CHANGf 

K = N 
100 K = K - 1 

IF(K,LE.O) RETURN 

I = . L ( K ) 

JFl I,LE#K» GO TO 120 

00 110 J = 1,N 

WAIT = X(J,K) 

X( J,K) = -X( J,I ) 
110 X(J,n = WAIT 
120 J = ^'(K) 

IF(J.LE.K) GO TO 100 

00 130 I = liN 

WA I T = X ( K , n 

x(K, n = -xt j,n . _j 

130 X(J,I) = WAIT 

GO TO 100 
600 lERR =1 

RETURN 

END 
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0001 SUBRGUriNfE C MPR S S ( NC , NEW , NC tN* N , \'V ) 

0002 DIMENSION D AT I X ( 1 00 , 10 0 > , POL < 2 5 , 1 0 0 ) , N V ( } - ) , '■'c ( 1 DO ) 

0003 COMMCN/OEP/OATI X»FMP,POL 

0004 DO 100 IC = 1 ,NC,3 

0005 [NC = NV(IC> 

0006 ! B = NV ( IC +i > 
J007 IFD = IMV'{ rC ■»-2> 
0008 IMV = NVdC +3! 
a009 I eD2 = NV( IC 

0010 DO 50 ! = 1,N0 

0011 SUM = 0.0 

0012 SUMP = 0.0 

001 3 00 AO J = IB, leo 

OOl-!* IF(MN.EQ.O.OR.MN.LT. n GO TO 40 

0015 SUMO = SUMD + POL (I, J) 

0016 W SUM = SUM + DATIX(I»J) 

0017 IF(MN.£0.0,OR.MN, LT. r ) GO TP 50 

0018 POL ( I , I NC ) = SUMP 

0019 50 DATIX(I,INC» = SUM 

0020 IF( IED+l.EQ.ieD2> GO TO 100 

0021 lA = lEO + 1 

0022 1 7. = IED2 - 1 

0023 IF{ IZ.LT.O) GO TO 100 

0024 DO 80 I = 1 ,N0 
002 5 K = INC + 1 

0026 00 RO J = I A, IZ 

0027 IF(MN,eO,0,aR.MN,LT. n GO TO TQ 

0028 PQL( I ,K) = POLd ,J) 

0029 70 DATrX(I,KJ = OATrX(I,J/ 

0030 80 K = K + 1 

0031 100 CONTINUE 

. 0032 IFdED, EO.NO) GO TO 25 

0033 DO 20 I = l,NC 

3034 K = INC+1 

0035 M = lED *\ 

0036 00 20 J = M,NO 

0037 IF(MN,E0.O,OR.MN, LT, 1) GO TO 19 

0038 . POL( I ,K) = POL( I ♦ J) 

0039 19 DATIX(I,K» = OATIX(I,J> 
004 0 2 0 K = K + I 

0041 25 IZE = NEW + 1 

0042 00 30 I =1 ,N0 

0043 00 30 J= IZE, NO 

0044 IF(MN,EQ.O,OR,MN,LT,I ) GO TO 30 

0045 POL ( I , J ) = 0. 0 

0046 30 DATI X{ I fJ) = 0. 0 

' 004 7 DO 2 00 IR = 1 ,NC,3 

0048 INR = NV(IR) 

, 0049 IB = NVCIR +1) 

0050 lED ,= NVdR +2» 

0051 IMV = NVdR +3) 

0052 IED2 = NVdR +4) 

0053 DO 170 J = 1,NEW 

0054 SUM = 0.0 

0055 DO 150 I = 18, 1 ED 

0056 150 SUM = SUM * OATlXd,J) 
0"')5 7 170 0ATI/X( I NRt J I = SUM 

ERI058 IF{ IED+I.EQ.IED2) GO TO 200 
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0059 lA = lED +1 

0060 II = 1602 -1 

0061 IF( IZ.LT.O) GO TC 200 

0062 DO 180 J = 1,NEW 

0063 K = INR > 1 
.0064 DO 180 I = IA,IZ 

0065 DATIX(K,J) = DATIXdtJ) 

0066 180 K = K +1 
.0067 200 CONTINUE 

006 8 IF{ I ED. EO.NC) GO 125 
J069 00 120 J = 1,NEW 

007 0 K = INR + 1 

0071 M = I ED + 1 

0072 00 120 I = M,Nn 

0073 DATIX(K.,J> = DATIXUtJ) 

0074 1 20 K = K+ I 

0075 125 !ZE = NEW +1 

0076 00 1 30 J = 1, NEW 

0077 00 130 I = IZE,NC 

0078 1 30 OATI X( I ,J) = 0. 0 

0079 00 3 00 IV = 1 ,NC ,3 

0080 INC = NV( rvi 

0081 IB = NV( IV +1) 

0082 lEO = NV( IV +2> 

0083 IMV = NVilV +3» 

0084 IE02= NV( I V +4) 

0085 SUM =OsO 

0086 DO 240 J = I8,IEn 

. 0087 240 SUM = SUM + EyP(J) 

0088 PMP{ INC ) = SUM 

0089 IF( IED+1.,EQ, 1ED2)- GO TO 300 
' 0090 lA = lED 1 

0091 IZ = IE02 - 1 

0092 IF(IZ«LT.O» GO TO 300 
009 3 K = INC + I 

0094 00 280 J = IA,IZ 

0095 EMP(K) = EMP(JJ 

0096 2 80 K = K ♦ 1 

0097 300 CONTINUE 

0098 IF( lEO. EO.NG) GO TO 225 

0099 K = INC + 1 

0100 M = lED + 1 

0101 00 220 J = M,NO 

0102 EMP(K) = EMP(J) 

0103 220 K = K + 1 

. 0104 225 I ZE = NEW + 1 

0105 00 330 J = IZE,NO 

0106 330 EMP( J) = 0.0 
- 0107 500 CONTINUE 

0108 RETURN 

0109 ENO 
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RLK DATA 



0 \ - 



- 7 : 7 



10/ 



0001 
0002 
0003 
0004 
0005 



BLOCK DATA 

UIV!E\SI0N OATA( 1 CO, J 00) ,PCL{ 2"^, 1)0) , f^MP M iO' 
COMMON/OEP/DATA ,EMP ,PCL 

OATA DATA/ 10 00 0*0. 0/ ♦ dOL/250 0*0. 0/ , f-^-^P/ 00'^ 1.1/ 
END 
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